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Abstract 

We show that the planar Chern-Simons (CS) theory on can be described by its 
dimensionahy reduced model. This description of CS theory can be regarded as a 
novel large- reduction for gauge theories on S^. We find that if one expands the 
reduced model around a particular background consisting of multiple fuzzy spheres, 
the reduced model becomes equivalent to CS theory on in the planar limit. In 
fact, we show that the free energy and the vacuum expectation value of unknot 
Wilson loop in CS theory are reproduced by the reduced model in the large-iV 
limit. 
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1 Introduction 

The large-A^ reduction [1] asserts that large-A^ gauge theories on flat space-times are 
equivalent to matrix models (reduced models) that are obtained by dimensional reduction 
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to lower dimensions (for further developments, see [SHS]). It is not only conceptually 
relevant because it realizes emergent space-time in matrix models, but also practically 
relevant because it can give a non-perturbative formulation of planar gauge theories as 
an alternative to lattice gauge theories. It is well-known that the equivalence does not 
hold naively due to the U{lY symmetry breaking [2]. Some remedy is needed for the 
equivalence to indeed hold (for recent studies, see [714TT]). However, as far as we know, 
there has been no remedy that manifestly preserves supersymmetry in gauge theories. 
Recently, the large-iV reduction was generalized to theories on in p23 (for earlier 
discussions, see [T314T5] ). Here is constructed as a nontrivial U{1) fiber bundle over S"^ 
by expanding matrix models around a particular background corresponding to a sequence 
of fuzzy spheres. This novel large- reduction is important from both of the above two 
viewpoints. It would give hints to the problem of describing curved space-time (TB] in the 
matrix models p!7HT9] . Unlike on flat space-times, the equivalence would hold on curved 
space-times without any remedy since the theories are massive due to the curvature of 
S^. The novel large- reduction can, therefore, give a non-perturbative regularization of 
gauge theories that respects both supersymmetry and the gauge symmetrjt. 

The novel large- iV reduction has been studied mainly for two cases so far. One is A/" = 4 
SYM on Rx p^[2T] - [23] . The reduced model of this theory takes the form of the plane 
wave matrix model [23]. The large- iV reduction provides a non-perturbative formulation 
of A/" = 4 SYM, which respects sixteen supersymmetries and the gauge symmetry so that 
no fine tuning of the parameters is required. Thus the formulation gives a feasible way 
to analyze the strongly coupled regime of TV = 4 SYM, for instance, by putting it on a 
computer through the methods proposed in [251427] . and therefore enables us to perform 
new nontrivial tests of the AdS/CFT correspondence [25H5D] . 

The other is Chern-Simons (CS) theory on S^, which has been exactly solved [31] . 
In this case, the equivalence between the reduced model and the original theory can 
be verified explicitly, which was briefly reported in [32] • In this paper, we present the 
verification of the equivalence in detail. We also discuss the large- iV reduction for CS 
theory on S^/Zg. Our result is an explicit demonstration of the novel large- reduction. 
CS theory on is a topological field theory associated with the knot theory and therefore 



^For recent developments in lattice approach to supersymmetric gauge theories, see [20] and references 
therein. 
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interesting in its own right. Our formulation gives a new regularization method for CS 
theory on S^. CS theory on is also interpreted as open topological A strings on T*S^ 
in the presence of D-branes wrapping S^. Topological strings capture some aspects of 
more realistic string, so that our study is expected to give some insights into the matrix 
models as a non-perturbative formulation of superstring. 

More recently, it was shown in [33j that the large- reduction holds on general group 
manifolds. In the case of ~ SU{2), the background taken in the large- A^ reduction 
in [33] is different from that in ^12j . It is desirable to explicitly demonstrate this different 
type of the large- A^ reduction, for instance, in the case of CS theory on S^. We leave this 
for a future study. 

This paper is organized as follows. In section 2, we briefly review CS theory on S^. In 
section 3, we summarize the relationships among CS theory on S^, 2d YM on 5*^ and the 
matrix model which were obtained in [3ll|35]. In section 4, we argue how CS theory on 

is realized in the matrix model as a novel large- A^ reduction. In section 5, we show the 
equivalence between the theory around a particular background of the matrix model and 
CS theory on S^. In section 6, by using the Monte Carlo simulation, we give a check of 
this equivalence. Section 7 is devoted to conclusion and discussion. In appendices, some 
details are gathered. 

2 Chern-Simons theory on 

In this section, we review some known facts about CS theory on needed for this 
paper [36l|37]. We consider CS theory on with the gauge group U{N), whose action is 
given by 



defines a topological invariant of the manifold 5''^, which also depends on a choice of fram- 
ing. Given an oriented knot /C in S*^, one can consider the Wilson loop in an irreducible 




(2.1) 



where k must be an integer. The partition function 




(2.2) 
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representation R of U (N) 

WRi]C) = TTR{Pexp (p A] . (2.3) 



The expectation value of the Wilson loop 

{Wr(}C))cs = ^ [vA Wr{}C) e''^^ (2.4) 
^cs J 

defines a topological invariant of K, depending on a choice of framing as well. In this paper, 
we are concerned with the Wilson loop for the unknot, which we denote by lyij(unknot), 
and mainly consider the Wilson loop for the unknot in the fundamental representation, 
which we denote by PVn (unknot). 

is obtained by gluing two solid tori along their boundaries with the transformation 

K = T"'ST'', (2.5) 

where 5* and T are generators of SL{2, Z), and m and n are arbitrary integers. T has the 
following representation in the Hilbert space 'H(T^), which is obtained by doing canonical 
quantization of CS theory on R x T^: 

{R\T\R') = Srr'Tr, (2.6) 



Tr = exp 



A-P + -A^ - —k{N^ - 1) 



(2.7) 



_k + N \ ^ 2 24 
where \R) denotes a state associated to the highest weight A, and p is the Weyl vector of 
SU{N). The partition function is given by 



■'CS 



m\o), (2.1 



where stands for the trivial representation, while the expectation value of the Wilson 
loop for the unknot is given by 

(H^^(unknot))c5 = (2.9) 

The canonical framing corresponds to m = = 0. In the following, we first present 
the exact results for the partition function and the expectation value of the Wilson loop 
in the canonical framing. The partition function fl2.2p was computed exactly: 

z-r-' = {o\s\o) = j^^^ E -p {-ytn' ■ ^^^0 ' ('-'^^ 



where the sum over w is a sum over the elements of the Weyl group W of U{N), and e{w) 
is the signature of w. The expectation value of the Wilson loop for the unknot was also 
computed exactly: 



2iTi 



One can obtain different expressions for (]2.10p and ( ]2.1ip by using Weyl's character 
formula 



and Weyl's denominator formula 

a ■ u 

I smh 

where a > are positive roots, and M/j is the set of weights associated to the represen- 



e{w)e"^P^-'' = JJ 2 sinh (2.13) 



tation R. By applying (12.131) to (I2.10p . one obtains 



27ri a ■ p 



Z-=-=o = i-— rr 2 sinh , (2.14) 



while by applying fl2:T2ll and f l27[3l) to flZTTl) 



For the fundamental representation, it reduces to 



TTiN 

(iyn(unknot))™r=° = ^^^^^S^- (2-16) 

smh^ 

Furthermore, by using (I2.13p . one can show that (12.101) is equivalent to an integral 
over N variables which takes the form analogous to a partition function of a matrix 
model [36l[38l|39]: 

jj_^e ...^^ n 

1=1 i<j 



sinh 
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where we have introduced 

- (2.18) 

which is identified with the string couphng in the context of topological strings. We 
call the statistical model defined by f l2.17p the Chern-Simons matrix model (CS matrix 
model). Correspondingly, one can find a relation 

(Tyfl(unknot))J^r=° = e-'^^^P-^+-^^'\chR{P))csM = e-S'^P-^+-^^'^ / ^ e'^'' 

(2.19) 

where 

<...w^^^ /nf (••■)--^'n(--^)^. 

Finally we consider the framing with m = n = 1, which appears in the direct evaluation 
of the partition function of CS theory on in [40J. One can see from fl2.8l) that the 
partition function in this framing is related to the one in the canonical framing as 

<j-m=n=l rji2 •7m=n=0 

^C5 — ^CS 

e-f(^^-i)+ff^(^^-i) r 2m a-p 

— 112 smh ' 



(k + N)^/^ Ai yk + N 2 

1=1 i<j ^ / 

One can also see from f l2.9p that the expectation value of the Wilson loop is related to 
the one in the canonical framing as 



(iyR(unknot))^r=' = ^(W^R(unknot))^r=° = {cMP))csm = ( Y1 ""'A 



CSM 

(2.22) 



In particular, because A ■ p + = ^'^^"^ for the fundamental representation, one obtains 

(iyn(unknot))-r=^ = (2-23) 
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Similarly, the perturbative part of the free energy is evaluated in this framing as 

N CO , -2 2 \ 

Tcs = %NiN^ - 1) + - E 1 + • (2-24) 

j=l n=l ^ ^ 

In the following sections, we will show that the reduced model of CS theory on 
reproduces the results for the original theory in the framing corresponding m = n = 1. 



3 CS theory on S^, 2d YM on 5*2 and a matrix model 

In this section, we review part of the results in [3ll[35], which we use in this paper. In 
section 3.1, we dimensionally reduce CS theory on 5*^ to 2d YM on 5*^ and a matrix 
model. In section 3.2, we describe a classical equivalence between the theory around each 
multiple monopole background of 2d YM on S"^ and the theory around a certain multiple 
fuzzy sphere background of the matrix model. In section 3.3, we exactly perform the 
integration in the matrix model. In section 3.4, using the result in section 3.3, we show 
that the equivalence in section 3.2 also holds at quantum level. 

3.1 Dimensional reductions 

In order to dimensionally reduce CS theory on S^, we expand the gauge field in terms of 
the right-invariant 1-form defined in (1A.3I) as 

A = iXiE\ (3.1) 

We choose an 5*0 (4)-iso metric metric for 5''^ flA.6|) fixing the radius of to 2//i. Then, 
we rewrite (12. ip as 



Scs = -^J JJ^^' (^tfie'''X,CjX, + /iXf + je'^'X,X,X,^ , (3.2) 



where £j is the Killing vector dual to defined in flA.7p . and we have used the Maurer- 
Cartan equation ( lA.SI) . 

As we explain in appendix |X1 we regard as a U{1) bundle over S"^. Here the fiber 
direction is parametrized by y. By dropping the y— derivatives, we obtain a gauge theory 
on 

Sbf = TTY- f ^Tr ( if^e'^'^X.L.X, + ^iXf + '^e^^'X,XjXk) , (3.3) 



where Qym = "ffc^' radius of S'^ is and Li are the angular momentum operators 
on S'^ given in (1A.13P with q = 0. qym will be identified with the coupling constant of 2d 
YM on 5"^ below. First, we see that (13. 3 p is the BF theory with a mass term on S'^. We 
expand as [T5|l^ 

X = ix { xer + aee^ r—^a^eel, (3.4) 

\ sva.6 J 

where Cr = (sin cos(y9, sin sincyj, cos6') and = = are the gauge field 

on S*^, and x is a scalar field on S"^ which comes from the component of the fiber direction 
of the gauge field on . Indeed, we can rewrite (13. 3p as 

,4 



Sbf = ^I ^Tr [xe^^f,^ - x') , (3.5) 



where f^u = d^a^ — dua^ + i[a^,ai,] is the field strength. This is the BF theory with a 
mass term: the first term is the BF term and the second term is a mass term. Next, we 
integrate x out in (13. 5p to obtain 2d YM on S"^: 

Furthermore, by dropping all the derivatives in (13. 3p and rescale Xj as Xj — > /iXj, we 
obtain a three-matrix model: 

Sm = -^Tr (xf + '-e'^''X,[X^,X,]] , (3.7) 

where Xj g^ = the sense of the Dijkgraaf-Vafa theory [42j, this matrix 

model is regarded as a mass deformed superpotential of A/" = 4 SYM, which gives the 
so-called N = 1* theory. 

3.2 Classical equivalence between 2d YM on S'^ and the matrix 
model 

The matrix model (13.70 with the matrix size M x M possesses the following classical 
solution, 

X, = L, = 04^'=1®U (3.8) 



^ While k in (|2.ip must be an integer, such a restriction is not imposed on k in p.3p . 
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where L^^' are the spin js representation of the SU{2) generators obeying [L^^', L^^'''] = 
ieijkL^l'\ and the relation + = M is satisfied, s runs over some integers, but 

its range is not specified here. It will be specified later as —A/2 < s < A/2 with A an 
even positive integer. 

We put 2js + 1 = Vt + Hs with f2 and ns integers and take the limit in which 

f2 — oo with — = — ^ — - = fixed, (3.9) 

where A = Att/^"^ is the area of S"^. Then, it was shown classically in [31] that the theory 
around (13.81) is equivalent to the theory around the following classical solution of (13. 3p 
with the gauge group U{K), 



fiL, + X, = /idiag(- ■ ■ , ■ ■ ■ , ■ ■ ■ ), (3.10) 



where Qs = ns/2, and Ng = K is satisfied. L^'^''' are the angular momentum operators 
in the presence of a monopole with the monopole charge g^, which are given in ( 1A.13I) . 
This theory can also be viewed as the theory around the following classical solution of 

(ES]), 




(3.11) 

where the upper sign is taken in the region < < vr and the lower sign in the region 
< < TT, and ag and represent the monopole configuration. Equivalently, this 
theory can be viewed as the theory around the multiple monopole background with the 
background for the gauge field given in (13. lip of 2d YM on 5*^ (13. 6p . As reviewed in the 
following two subsections, the above equivalence is extended to quantum level 

3.3 Exact integration of the matrix model 

First, we rewrite the matrix model (13.70 in terms of the following matrices; 

Z = Xi+iX2, 
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(3.12) 



where Z is an M x M complex matrix and $ is an M x M hermitian matrix. The partition 
function is defined by 

^ra = UT]^.^^^ [ dZdZU^e'^^, (3.13) 



vol(f/(M)) 

Sm = -^tr Z^] + (1 - it)ZZ^ + (1 - le)^^] . (3.14) 

9m 

Here, in order to make the path integral well-defined and converge, we introduce i in 

the exponential and —ie in the action. Taking the gauge in which $ is diagonal as 
$ = diag(0i, 02, ■ ■ ■ , 0m) and integrating Z and Z"^ out, then we obtain 

- (2^0^ y + ^"^ ' (3-15) 

where Hj^jl'/'i ~ 'Pj) comes from the diagonalization of $, which is the square of the 
Vandermonde determinant. (27r)^^ is the volume of f/(l)^^ and M! is the volume of the 
Weyl group of ?7(M). 

Next let us concentrate on the factor in the integrand in (I3.15p . 

■ (3-16) 

-L -1 0. - . + 1 _ 

For this factor, we use the following formula, 

^ P.V.- + m5{x), (3.17) 



X — le X 
where "P.l/." denotes the principal value. Then we obtain 



TT Vj 

\}: 0i - 0j + 1 - ie 



n 



- (j)j){(f)j - (f)i) { iT^\^f^ ^ , i\ , 



(3.18) 

It is easily seen that f l3.18p is written in the sum of terms containing 



111 
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2 y 



2j2 ^^2j2 + l 



f ) ' m"" - - - ^r^ - 1) ■ ■ ■ m;^ - 0^2^+1 - 1), (3.19) 

where K can take the value 1, 2, ■ ■ ■ , M and we have reordered and relabeled the eigen- 
values of $ as 

$ = diag(0«, ■ ■ ■ , ■ ■ ■ , ■ ■ ■ , <f>f\ ■ ■ ■ , (3.20) 

The size of the matrices, M, satisfies M = Ylf=ii'^ji + '/'l^^ represents the i-th 
component of the J-th block. Thus (13.181) is decomposed into terms, each of which is 
characterized by an M-dimensional reducible representation of the SU{2) algebra con- 
sisting of K blocks (irreducible representations), with a U{1) degree of freedom in each 
block. We label the reducible representations by "r"|f| and denote the U{1) part in the 
J-th block by aj, putting aj = (t>2ji+i + ji- Then, we obtain 



TT TT TT — — — t 1 <^'^'' ■---^=^--™/=-^i 

I^J mj=-jj nij=-jj 



ai + mj-aj- rrij --^(i-i.)Ef=iEi^=-,>K+"^i)' 



(3.21) 



where 



TV"^ = TT \ , (3.22) 



(tl of blocks with the same length)! 
After some calculations, we find that f l3.2ip results in 

{aj - aj)2 - {jj - jjf --|r(l-i^)Ef=i{2i7+l)a2 



K 



X />.K y n rfa, n + ^.23) 

Thus the partition function of the matrix model is decomposed into sectors characterized 
by the (reducible) representation of the SU(2) algebra. Note that P.V. prevents the 
integration over the U{1) part aj from mixing the sectors. 



Namely, "r" specify the number of irreducible representations of SU{2), K, and the dimensions of 
the irreducible representations, 2ji + 1 (/ = 1, • • • K), which satisfy the relation M — ^f^ii^ji + 1). 
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3.4 2d YM on 5^ from the matrix model 

In this subsection, we show that the partition function of YM on S"^ is obtained from 
that of the matrix model [52]. As we will see, the sector which consists of K irreducible 
representations in the matrix model corresponds to SU{K) 2d YM on S"^. In the K block 
sector, we group the irreducible representations with the same dimension and label the 
groups by s, such that js 7^ jt for s ^ t. We denote the multiplicity of the s-th group 
by Ng. Then, the relations K = Xls^s M = '^gNs{2js + 1) are satisfied. We also 
denote the i-th U{1) part in the s-th group by a^j, where i = 1, - ■ ■ , Ng. Note that in this 
notation TV; in (Km equals 1/(11. ^s!)- 

Since the partition function of the matrix model fl3.23p is completely decomposed into 
sectors without overlap, we can extract the sectors consisting of K blocks. In the i^'-block 
sector, due to the factor exp(— :A-(1— ze) J2f=^ i/(j/ + l)(2j/ + l)), configurations of almost 
equal size blocks are dominant in the limit in which M/K = Q ^ 00, g^/^ = fixed. So, 
we consider configurations around the dominated configuration by setting 

2js + l = n + Us, (3.24) 

where \ns\ ^ Q. Substituting this into the i^'-block sector in (13.231) and taking the limit 
( 13. 9p . we obtain 

2^/<-biock =c ^ / n n ^^"^ n n n ] ~ '^^^^^ ^ 4*^^^* ~ i 

{ns,Ns} s i=l s<t i=l j=l ^ ' 

i y-'Vs ( .2, 

X e ^YM^ V / . (3.25) 

Here we have absorbed overall irrelevant constants and divergences into a renormalized 
constant C. We have also replaced the denominator in the integrand in (I3.23P by constant 
(fi^) and taken the integral regions of asi to (—00, c)o), since the exponential factor in (I3.23P 
oscillates rapidly around |asj| > Q in the limit where e — )■ and f2 — t- 00, so that the 
integral dominates around \asi\ ^ ^l. 

Finally, by rescaling asi by ysi = 2asi, and making an analytical continuation Qym ~^ 
—iQYM^ we obtain 

„ Ns Ns Nt 

{ns,Ns} s i=l s<t 1=1 j=l 

12 



X e 



(3.26) 

where irrelevant constants are again absorbed into a constant C . Zk exactly agrees with 
the partition function of U{K) 2d YM on S"^ ptii7]H In YM on 5"^ are identified with 
the monopole charges since the dependent term in the exponential coincides with the 
contribution of the classical solution (13. lip to the action of 2d YM. This fact is consistent 
with the classical equivalence reviewed in section 3.2. 

4 CS theory on from the matrix model 

In this section, we provide a prescription for realization of CS theory on in 2d YM on 
S"^ and in the matrix model. We also generalize the prescription to CS theory on /Zq. 
In section 5, we prove the large-A^ equivalence between CS theory and the matrix model 
based on this prescription. 

CS theory on S'^ is realized from 2d YM on 5*^ through the large- reduction for 
the case of compact space proposed in [12] as follows. Note that the Kaluza-Klein (KK) 
momenta along the fiber (5*^) direction in are regarded as the monopole charges on S'^. 
Hence, if one considers the particular monopole solution in 2d YM on S*^ such that the 
spectrum of the monopole charges reproduces the spectrum of the KK momenta along , 
the theory on S*^ expanded around such monopole solution is equivalent to CS theory on 
S^. More precisely, in order to realize CS theory on S^, we first put = s and Ng = N 
in fl3.10p (or in (13. lip ) and then make s run from —A/2 to A/2, where A is a positive 
even integer. Finally we take the following limit, 

A^oo, AT^oo, -^k^ = = fixed. (4.1) 

27r A; + A^ ^ ' 

Then, we see that 2d YM on S"^ expanded around the monopole background is equivalent 
to U{N) CS theory on S"^ in the planar limit. Here, we can identify the difference of the 
monopole charges ng — nt = s — t with the momenta n & Z along the fiber direction and 
A is interpreted as the momentum cutoff. Introducing N and taking the limit N ^ oo 



fact, only the SU{K) part agrees, and the U{1) part does not (see ref. [35]). However, this 
difference does not matter in the fohowing arguments. 
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are needed to suppress the non-planar contribution [12]. As for the relation between the 
coupling constants of the two theories, we naively expect that — ^^2-^^ ~ T ~ fixed. 
However, the last relation in (14. ip implies that we need some renormalization. In sections 
5 and 6, we show that the above prescription including the renormalization of the coupling 
constant is indeed valid. 

Furthermore, by combining this equivalence with the relationship between 2d YM on 
S"^ and the matrix model reviewed in section 3.4, we obtain the following statement: if 
one takes in (13. 8p 

A A 

< s < -, 

2 - - 2' 

2j, + i = n + s, 

Ns = N (4.2) 

and takes the limit in which 

i7 — 7- cxD, A — )■ cxD, f2 y oo, N ^ oo, 

2 

^ = ' = fixed, (4.3) 

the theory around (13.81) of the matrix model (13. 7p is equivalent to U{N) CS theory on 
5*^ in the planar limit. In sections 5 and 6, we confirm this statement. The particular 
background (14.20 and the limit (14. 3 p is the same as the ones adopted in realizing Af = 4 
SYM on Rx in PWMM [12]. 1^ and A play the role of the ultraviolet cutoffs for the 
angular momenta on 5*^ and in the fiber direction, respectively. 

The results in section 3.4 and the above prescriptions lead us to consider the following 
statistical model, 

A/2 N N 

^= / n n^^- n ui^y^^^-y^^^'+^'-tne-^^^--/^^^-^^-. (4.4) 

s=-A/2 i=l -A/2<s<t<A/2 j,i=l 

This model is obtained from (I3.26P by making an analytic continuation as y^j Wsi, 
—7- —g"^ and putting Us = s, Ng = N {s = —A/2, ■ ■ ■ , A/2) . It follows from the above 
arguments that this model reproduces the planar limit of CS theory on in the limit in 
which 

A oo, oo (4.5) 
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with 




2TiiN 



QsN = fixed. 



(4.6) 



k + N 



In the next section, we will verify the equivalence of CS theory on 5*^ with (14. 4p which is 
equivalent to the reduced matrix model expanded around the particular background. 

We can easily extend the above argument to CS theory on S^/Zq. We consider the 
case in which Zg acts on the circle along the fiber direction in 5*^. If we use the notation 
in appendix A, the periodicity for S^/Zg is expressed as, 



so that the radius of the fiber direction is 1/g of that for q = 1 and the momenta along 
the fiber direction are discretized as 0, ±g, ±2g, ■ ■ ■ . Then, we see that if one imposes 
(14. 2 p with the second equation replaced with 2js + 1 = Q + qs and takes the limit (14. 3p . 
the theory around (13. 8p of the matrix model (13. 7p is equivalent to U{N) CS theory on 
S^/Zq in the planar limit. Correspondingly, the statistical model for CS theory on S^/Zq 
is given by replacing s, t in (14. 4p with qs, qt. 

5 Proof of equivalence 

In this section, we prove our statement that the model (14. 4 p reproduces the planar limit of 
CS theory on in the limit shown in (14. 5 p and (14. 6p . This equivalence can be understood 
from a correspondence between the Feynman diagrams in the two theories. The Feynman 
diagrams in (14. 4p have a one-to-one correspondence to those in (I2.17P which was obtained 
by rewriting the partition function of CS theory on S^. Furthermore, if one takes the 
limit given by (14. 5 p and (14. 6p . each Feynman diagram in (14. 4p takes the same value as the 
corresponding diagram in (I2.17p . We consider the free energy and the vacuum expectation 
value of unknot Wilson loop operator in the model (14. 4p to see the equivalence. We show 
that they coincide with the results in CS theory on completely in the planar limit. 

5.1 Feynman rule for our matrix model 

Let us first consider the Feynman rule of (14. 4p . For this purpose, we rewrite (14. 4p into a 
manifestly U{N) invariant form. We will see that the result is given by a multi-matrix 



{9, if, ^ {e,ip + 27r, + 27r) ~ {6, ^ + 47r/g) 



(4.7) 
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model with double trace interaction. This form allows us to compare the two theories 
explicitly. 

We first factorize the measure in (14. 4p as follows, 

niKfc. - y^f + - = E(-' - ir'Wy- - i-^)'!!!! (i + 

s<t i,j s<t s i<j s<t i,j ^ \ / / 

(5.1) 

The first factor on the right hand side is just a constant which we will omit in the following 
and the second factor gives the Vandermonde determinant for each s. We next introduce 
the 't Hooft coupling in this model as, 

A = 27T'^g^N, (5.2) 

and redefine the fields as (J) si = 2iTysi/\^. Then, the third factor in (15. ip can be written 

as 

s<t i,j ^ ^ ' 

-P<l2 2l^Li°g^i + i^^7T7)^ 

sj^t i,j ^ ^ ' 

oo / , \ fc 2fc 



exp 



Finally, defining U{N) covariant matrices 0s such that (psi is i-th eigenvalue of 0s for each 
s, we find that (14.41) can be written as a multi- matrix model. 



n dcj>s exp I -y J] tr02 - \/(0s) I , 



Z = y _[|(i0sexp|--^tr0^-y(0s)}>, (5.4) 
where V^(0s) is given by the double trace interaction, 

^■w.) - EE^^ (i^^)'E ( Z ) -cM-w-". (5.5) 

s^t k=l ^ / \ / jn=o \ / 

The Feynman rule for (15. 4p is given as follows. The tree level propagator is 

{(psij(l)tki) = j^^ii^jkSst- (5.6) 
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10 



Figure 1: The propagator (left) and the interaction vertex (right) in the model fl4.4p . The 
right figure shows the vertex which includes tr0^tr0j . The number 10 represents the total 
number of 0s and (pt included in the vertex. 
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Figure 2: Typical connected Feynman diagrams which contribute to the free energy in 
the model (15^ . 



In terms of the standard double line notation, this is represented by double lines with 
indeces s and t as shown in Fig{Tl The double trace interaction vertices are represented 
in terms of dashed lines which connect two traces. For example, FigH] (right) shows the 
term in fl5.5p which includes tKp^trcp'l . We connect two vertices coming from tr0^ and tKpf 
in terms of a dashed line to represent a double trace interaction vertex. The weights of 
the vertices can be easily read off from (15.51) . The tr0^tr0j vertex gives a constant factor. 



1 



) a + b / 
' (-i)M " + ' 



(5.7) 



As in the standard perturbation theory, physical observables such as the free energy 
are computed by summing all the connected diagrams. Figj2] shows typical examples of 
connected diagrams which contribute to the free energy. More precisely, the connected 
diagrams in this theory are those which are connected by dashed lines or by double lines. 
One can show that only such diagrams indeed contribute to the calculation of observables. 

Since we are interested in the planar limit, we consider what type of diagrams is 
dominant in this limit. It is easy to see that the leading order in the 1/A^ expansion is 
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given by the diagrams that satisfy two conditions. One is that they are planar diagrams 
with respect to the double lines in the ordinary sense. The other is that they are divided by 
two parts by cutting any dashed line. We call the diagrams satisfying the latter condition 
'tree' diagrams in this paper. We can see that the diagrams in which the dashed lines 
form any loop (i.e. non-'tree' diagrams) give subleading contribution. In Figj2l the left 
diagram contributes to the free energy in the large- limit while the right diagram is 
suppressed since the dashed lines form a loop. 

5.2 Feynman rule for Chern-Simons matrix model 

Next, we construct a Feynman rule of CS matrix model f l2.17p . We first rewrite (12.171) 
into a manifestly U {N) invariant form as we did for (14. 4 p in the previous subsection. The 
identification for the coupling constants shown in (14. 6 p allows us to put 

A = g,N, (5.8) 

where A is the 't Hooft coupling in (15. 4p which is introduced in (15. 2p . If we redefine the 
fields as (pi = /3i/\/X. we can rewrite the measure in (I2.17P as follows, 

Tfr^ sinh^ 



n (2 sinh ^) ~ n('^^ - '^^•)' E ■ v5 

i<j ^ ^ i<j L i,j 2 



00 



n(0.-0,)'exp 5^5^1og(l 

i<j L i,j n=l 



n(0.-0,)^exp 5^5^ 

i<j L i,j k=l 




-cm 



\ 2k— m 



(5.9) 

where ~ represents that we have omitted an irrelevant constant factor and C(z) is the 
zeta function. We have used the formula, 

sinhvra; tt -.^n 

— =n(i+;^). (5-10) 

n=l ^ ' 

to obtain the second line in (15. 9p . If we introduce an U{N) covariant matrix such that 
its 2-th eigenvalue is given by (pi, we obtain the following partition function for CS matrix 
model on 

ZcsM ~ j d^exp |-ytr02 - ^(0)| , (5.11) 
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where the potential V{(j)) is given by the double trace interaction, 

f (^) -t^-iri^)'t(Z) trrtr(-^)--. (5.12) 

k=l ^ ^ m=0 \ / 

We construct the Feynman rule of ( 15. lip as follows. The tree level propagator is given 

by 

{(pij(pki) = jj:Sii6jk, (5.13) 

which is expressed by double lines as usual in the standard double line notation. The 
double trace interaction vertices are represented in terms of dashed lines which connect 
two traces as in the case of (15.41) . They are represented in the same way as in Fig JT] except 
that we do not need the indeces s, t here. Fig{T] (right) without s, t corresponds to the 
term in (15.121) which includes tr0^tr0^. The weights of the vertices can be easily read off 
from (15.121) . The tr^'^tr^'' vertex gives a constant factor, 

4C(a + 6) f-X\'^ , ( « + & 



a + h V 47r 



ir • (5.14) 



As in the case of (15. 4p . connected diagrams in this theory are those which are connected 
by dashed lines or by double lines. In the N oo limit, the leading order in the 1/A^ 
expansion is given by the diagrams which are planar in the ordinary sense and 'tree' with 
respect to dashed lines. For instance, the right diagram in Fig|2]in the case of (15. lip (i.e. 
the diagram without the indeces s,t,u, - ■ ■ .) is suppressed in the large- limit since the 
dashed lines form a loop. Note that the Feynman diagrams in (15. lip have the one-to-one 
correspondence to those in (15. 4p . 

5.3 Equivalence between two theories 

In this subsection, we show that (15. 4p and (15. lip are equivalent in the limit shown in (14. 5 p 
and (14. 6p . We can prove this equivalence by showing that, in this limit, each Feynman 
diagram in (15. 4p takes the same value as the corresponding diagram in (15. lip . Note 
that we have already imposed the identification of the coupling constants shown in (14. 6 p 
through and (ISlSj) . 
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5.3.1 Free energy 

In order to prove the equivalence, we first show that 



^ (5.15) 



A^2(A + i) ^2 

holds in the limit shown in (14. 5p and (14. 6p . where T and Tm are the free energies of 
(15 ■4p and (15. lip , respectively. As an example of the Feynman diagrams contributing to 
the free energy, let us consider the left diagram in Fig|2]for the both theories. Since the 
only deference between (15.70 and (15. lip is 2(s~tY+^ Cl*^ + it is easy to see that the 
diagram in (15. 4p divided by A^2(A + 1) becomes equal to the corresponding diagram in 
(imj) divided by if 

^ A/2 A/2 A/2 A/2 A/2 ^ 

24(A + 1) is^tfU-uY^{u-v)Hu-wY 

{tj^s) {u^t) {v^u) {w^u) 

= C(8)C(10)C(4)2 (5.16) 

holds. In fact, this identity holds in the A — > oo limit. We show the proof of (I5.16P 
in appendix |Bl Thus, we see that the two diagrams with the appropriate normalization 
shown in (I5.15P take the same value in the A — )■ oo limit. 

We can show more general formula which includes (15.161) as a special case. Let us 
consider generic 'tree' planar diagram in (15. 4p shown in FigJS] For this diagram, we can 
show the following identity in the large-A limitj^l, 

1 1 " 

2"(A + 1) Z (s,_,2)2fc.(,^_,3)2.....(,^_,^^^)2.„ = n (5.17) 

where {sa} represents the set of indeces (si, ■ ■ ■ ,s„+i) and each Sa runs from —A/2 to 
A/2. Note that the singular points such as si = S2 are not included in the summations as 
in the simpler case of (I5.16P . The proof of (I5.17P is shown in appendix |Bl (I5.17P implies 
that the value of the diagram in Figj3] divided by N'^{A + 1) is equal to the same diagram 
in (15. lip divided by A^^. We find, therefore, that (I5.15P holds to all orders in A in the 
limit fOjl . 



^ Note that ()5.17p holds also for the diagrams which are 'tree' and non-planar since (|5.17p depends 
only on the structure of the dashed lines. 
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Figure 3: Generic 'tree' planar diagram in (15 ■4p . Each blob represents a planar diagram 



We can check f l5.15p also by performing an explicit perturbative calculation. The detail 
of the calculation is shown in appendix O Up to 0{X^), the result is given by, 

+ ^ + 



N^{A + 1) 12 288 1440A^2 
A3 



+ 



I 9072 487r6A (s - t)^{t - - s)^ 
while for CS matrix model, we can read off the counterpart from (12.241) as, 

= 12 " 12m + 288 + 96A^ + ^^^'^^ ^^'^^^ 

Thus, we see that if one takes the limit shown in ( 14. 5 p and ( 14. 6p . ( I5.18P indeed agrees 
with (1519|) . 

5.3.2 Unknot Wilson loop 

In terms of the formula (I5.17p . we can show that other physical observables also take the 
same values in the two theories. In the following, we show that the VEV of unknot Wilson 
loop in CS theory which is given by (I2.23P is also reproduced from the model (14. 4p . 
We introduce the Wilson loop operator in A/" = 1* matrix model proposed in 

W{>C) = ^Tr 



"1 ^.Mi,, , , 

(5.20) 



Pexp(^^/ijr X,ElMa))^^^da^ 



where z'^ are coordinates on 5*^, cr parameterize the knot /C and E\,^dz^'^ are the right- 
invariant 1-forms on defined in ( IA.4p . In [M], it is shown through Taylor's T-duality 
that this operator is reduced to the Wilson loop operator on 5*^ in the continuum limit. 
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In order to see the correspondence to the unknot Wilson loop in CS theory, we consider 
fl5.20p with an unknot contour. We take the simplest one, a great circle on S^, as such 
contour. The great circle is parameterized, for example, by 



(^,V^,V')= (|,0,47ra), (5.21) 



with a G (0, 1]. By substituting f l5.2ip into fl5.20p and setting /i = 1, we obtain, 

ly (unknot) = ^Tr (e^"^^») . (5.22) 

As in section 3.3, we take the gauge in which X3 = $ is diagonalized. Recall that after we 
integrate out the delta functions shown in (13.191) which appear in the partition function, 
the remaining eigenvalues of $ are represented by ?/si/2. Then, in (I5.22p . if we make the 
analytic continuation i/si — )■ —iysi and the field redefinition y^j = V^0si/27r, we obtain, 

(unknot) = § E ^"^^^^ = mTTT) ^ ^^-^^^ 

s,i ^ ' s 

Therefore, we find that the unknot Wilson loop operator in the fundamental representation 
is given by in ( lOll . 

On the other hand, in CS matrix model fl2.17p . the VEV of unknot Wilson loop is 
given by fl2.22p . Then, the equivalence for the unknot Wilson loop is expressed in terms 
of and (ps as 

_i_ ^(Tl,v^*.) = ^(lVeV^*>„, (5.24) 

where (■ ■ ■ ) and (■ ■ ■ )m denote the expectation values in (15. 4 p and (15. lip , respectively. 
We can show (I5.24p in all orders in perturbation theory. In perturbation theory, the VEVs 
in (I5.24P are calculated by expanding the exponentials. Then, (I5.24p is satisfied if 

^EW:) = (Tr0")M (5.25) 

s 

holds for n = 0, 1, 2, ■ ■ ■ . Because (15.171) holds for all the 'tree' planar Feynman diagrams 
contributing to the VEVs in (15.251) . we find that (I5.25P and hence (15.241) hold in all 
orders in A. Thus, we have shown that the VEV of unknot Wilson loop in CS theory is 
reproduced from the reduced model. 
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It is straightforward to check (15.241) exphcitly by calculating the VEV of f l5.23p in 
perturbation theory. The method of calculating the VEV is almost the same as in the 
case of the free energy. We show the detail of this calculation in appendix O In the 
large- A limit, the result is given by 

TO - ^44-2^4^2!^^ (^■^'') 

On the other hand, (12.231) is expanded as, 

(H'o(uuknot))c. = l + A + + + O(A'). (5.27) 

Comparing (I5.26P and (I5.27p . we can see that (I5.24p indeed holds in the — )■ oo limit. 

5.4 Extension to S'^/Zq 

We can show the similar equivalence for the case of CS theory on S'^ /Zk- Although this 
theory possesses many nontrivial vacua, we consider the trivial vacuum only in this paper. 

The vacua in CS theory on /Zq are characterized by the holonomy along the circle 
on which Zq acts. In |l8], it is shown that the partition function of CS theory on /Zq 
for each vacuum sector is reduced to the partition function of a matrix model similar to 
(15. lip . For example, the partition function for the trivial vacuum sector is given by (15. lip 
with A replaced by A' = QsN/ q. 

On the other hand, as mentioned in the last part in section HJ the statistical model for 
CS theory on S^/Zq is given by replacing s,t in (14. 4p with qs,qt. By following the same 
calculation as (15. ip and (15. 3p . we can see that such model is equivalent to (15. 4 p with A 
replaced by A" = 2-7^"^ g"^ N / q^ . 

Then, with the identification A' = A", we see that the two theories are equivalent in 
the limit (14. 5p . Thus, we find that the theory on /Zq is also reproduced from the matrix 
model. 

6 Monte Carlo simulation 

The large- A^ equivalence between (14. 4 p and CS theory on 5*^ can be understood also from 
the agreement of the eigenvalue densities in the two theories. In this section, by performing 
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a Monte Carlo simulation, we show that the eigenvalue density of Usi at the saddle point 
of fl4.4p coincides with that of /3j in 02.1 7p if |s| is sufficiently small compared to the cutoff 
A/2. The eigenvalue densities for with s near the cutoff have some deviation from 
the density of However, we will see that such cutoff effect vanishes in the continuum 
limit and does not contribute to values of physical observables. From this property of the 
eigenvalue density, we can also give another check of the equivalences for free energy and 
unknot Wilson loop operator which we have shown in the previous section. 
Let us introduce the eigenvalue density of /3j in CS theory on 5''^ as, 

/^(^) = ]^E^(^-/^^/(2^))' (6.1) 
1=1 



and the eigenvalue density of Usi in (14. 4 p as. 



1 ^ 

Y,6{x-y,;). (6.2) 



i=l 

We also define the difference between (16. ip and (16. 2p as, 

ApW(a;) =p(^)(x) -p(x). (6.3) 

By performing a Monte Carlo simulation of (14.40 . we can see that ApW(x) = for 
s ~ 0{1) while Ap'^'^^(x) 7^ if s is near the cutoff ±A/2. Figure H] shows the result 
of the simulation. This result shows that as one goes toward the midpoint s = 0, the 
distribution converges rapidly to the distribution of CS theory on . From this result, 
we conclude that the deviation of the eigenvalue density exists only for sufficiently large 
\s\. 

We have also verified through the Monte Carlo simulation that ^p^^\x) possesses the 
following property, 

1 f 

hm— — dx J2 |Ap(^)(^)l=0. (6.4) 

s=-A/2 

This fact is shown as follows. Let us consider the saddle point equation of (14.40 . 

9 r f 1 11 

X- dx' V p^'\x') \ — + — } = 0. (6.5) 

J yx — x'^iys—t) X — x' — i{s — t) ) 



9' 



t=-A/2 
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Figure 4: The eigenvalue distribution of ysi in our model is plotted for s = 0, 28, 32, 36, 40 
with A = 0.25, A = 80 and = 48. The dashed line represents the distribution p{x) in 
CS theory on 5*^. The error bars are omitted since they are smaller than the symbol size 

Because ^p^^\x) = for s ~ in order for (16. 5 p to be satisfied, Ap*^*^(a;) with s near 

the cutoff should behave as 

Then, we can expect that 

A/2 

dx J2 |ApW(a:)|~^, (6.7) 

s=-A/2 

in the limit A — oo so that (16. 4p is satisfied. We confirmed (16.70 through the Monte 
Carlo simulation. Figure [5] is the result of the simulation and it shows that the quantity 
on the left hand side in (16.70 is indeed nicely fitted by f{x) = a + 6^2£A with a ~ 0. 

If we assume (16. 4p or (16. 7p . we check that the equivalence of the free energy (15.150 
indeed holds. Free energy of ( 14.40 is written in terms of ( 16. 2p as. 




(6.8) 
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Figure 5: J dxJ2s=-A/2\^P^^^(^)\ plotted against logA/A. A = 0.25 and has 
already been extrapolated from A^ = 24,32,48 to A^ = oo by a linear function of l/N"^. 
The dotted line represents the extrapolation of the plots using /(x) = a + b^^^^ with a 
and h constants. The error bars are smaller than the symbol size. 

Then, by substituting f l6.3p . we obtain 



By using the formula flS.lOp . we can show that the first line on the right hand side of 
fl6.9p is equal to J-'m/N'^ in CS matrix model in the limit A — > oo. We can further show 
that the second and the third lines vanish in A — t- oo as follows. Let us introduce a 
sufficiently large but finite constant b such that the eigenvalue densities for all s fit in 
the interval [—b,b]. Such b exists because the attractive force in (16.50 coming from the 
Gaussian potential is stronger than the repulsive force for |x| ^ 1 so that p^^\x) has a 
compact support. Then, the first term in the second line is bounded as 



AA^2 





(6.10) 
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Because of (16 ■4p . we find that this term does not contribute to J^/{AN'^) in the A — t- cxo 
hmit. We can apply the same argument to all the terms in the second and the third lines 
in ( 16. 9p . Thus, we have given another check of ( IS.lSp . 

We can also show the equivalence for unknot Wilson loop operator. Since the VEV in 
( 14. 4 p can be written as 

(ly(unknot)) = = N{A + i) / dxJ2{P^^\^)y''^ (6.11) 

Using ( 16. 4p . we find that 

(^(unknot)) ^ ^ y dx{p{x))e^'''' = (VTn (unknot ))cs, (6.12) 

in the limit A — t- oo. Thus, we have given another derivation of the equivalence for unknot 
Wilson loop. 

7 Conclusion and discussion 

In this paper, we studied a matrix model that is obtained by the dimensional reduction 
of CS theory on to zero dimension. We found that expanded around the particular 
background corresponding to S^, it reproduces the original theory on 5''^ in the planar 
limit. We calculated the partition function and the VEV of the Wilson loop in the theory 
around the background of the matrix model and verified that they agree with those in the 
planar CS theory on S^. We checked these results also by performing the Monte Carlo 
simulation. We also extended this equivalence to the case of CS theory on S^/Zg. 

In this paper, we only consider the unknot Wilson loop. It is relevant to see whether 
the VEV of the Wilson loop for a real knot is reproduced in our formulation. We should 
verify that the VEV of the Wilson loops does not change against continuous deformation 
of the loops, such that it indeed represents a topological invariant. We should also extend 
our analysis to the case of the Wilson loop with general representation TZ of SU{N). 
It is also interesting to see whether we can use our formulation to achieve the large- iV 
reduction for other theories including the CS term such as the ABJM model |49j. 
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A and 

In this appendix, we summarize some useful facts about and S"^ (See also 

is viewed as the SU (2) group manifold. We parameterize an element of SU{2) in terms 

of the Euler angles as 

g = e-*^'^3/2g-i9a2/2g-iV^3/2^ (A.l) 

where 0<^^<7r, 0<v9< 27r, < ^ < An. The periodicity with respect to these angle 
variables is expressed as 

{9, if, ^) r^{9,ip + 27r, + 27r) ~ {9, ip,i; + An). (A.2) 

The isometry of is 50(4) = SU{2) x SU{2), and these two Sf/(2)'s act on g from left 
and right, respectively. We construct the right-invariant 1-forms, 

dgg-' = -tixE'ai/2, (A.3) 

where are explicitly given by 

E-^ = —{— smipd9 + sm9 cosipdip), 
A* 

E"^ = — (cos Lpd9 + sin 9 sin Lpdif)) , 
/i 

E^ = -{d^ + coi^9dip). (A.4) 
/i 



The radius of is 2/yU. Ei satisfy the Maurer-Cartan equation 

2 



dE' - ^e,.jkE^ AE^ = 0. (A.5) 
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The metric is constructed from as 



ds^ = E'E' = \ {de"^ + sin^ Odif'^ + {dip + cos Odif)'^) . (A.6) 
/i 



The Kilhng vector dual to i?* is given by 



A = —EfdM, (A.7) 
/i 



where M = Oyipyip and El are inverse of E\^. The Kilhng vector is explicitly expressed 
as 



£i = —i ( — sin (fdg — cot 6* cos ipd^ H — ^ — Trdih) 
V smfc' / 



sm(/9 



£2 = — M cos o^c^e — cot 6^ sin o^c^io H — : — -d^ 
* sm 



£3 = -id^. (A.8) 

From the Maurer-Cartan equation (lA.Sp . one can show that the Killing vector satisfies 
the SU{2) algebra [Ci, Cj] = ieijuCk- 

Next, let us regard as a U{\) bundle over S'^ = SU{2) /U{1). S"^ is parametrized 
by 9 and (p and covered with two local patches: the patch I defined by < < tt and 
the patch II defined by < < vr. In the following, the upper sign represents the patch I 
while the lower sign represents the patch II. The element of SU (2) in flA.ip is decomposed 
as 

g = L-h (A.9) 

with 

h = e-^('/'±^)'^3/2^ ^A.IO) 

L represents an element of 5*^, while h represents the fiber U{1). The fiber direction is 
parametrized hj y = ip ± ip. The zweibein of S'^ is given by the i = 1,2 components of 
the left-invariant 1-form, —iL^^dL = /ieVj/2, which takes the form 

= — (± sin (pdd + sin 6 cos cpdip) , 
29 



= —{— COS ipd9 ± sin 9 simpdip). (^-11) 
1^ 



This zweibein gives the standard metric of S"^ with the radius 1//^: 

1 



ds"^ = e'e' = -^ide'' + sin' Od^'). (A. 12) 



Making a replacement dy —iq in (lA.SP leads to the angular momentum operator in the 



presence of a monopole with magnetic charge q at the origin [50] : 

r(<?) ■( ■ A , A ^ 1 Tcos^ 

L\ = t[sm (fOe + cot 6/ cos o^a,. ) — q — coso?, 

sm 

t(i) ■( A , ■ ^ \ 1 Tcos^ . 

-^2 = COS (poe + cot usmipOu,) — q — smy?, 

sm 

if = -id^Tq. (A. 13) 

where q is quantized as q = 0,±|,±1,±|,---, because y is a periodic variable with the 
period Air. These operators act on the local sections on S'^ and satisfy the SU (2) algebra 
[L'f\L^j^^] = ieijkLf^ . When g = 0, these operators are reduced to the ordinary angular 
momentum operators on (or Li = Lf\ The SU{2) acting on g from left survives 
as the isometry of S"^. 



B Proof of equation (15.171) 



General 'tree' planar diagram shown in FigJ3] consists of several parts (blobs) which are 
connected by the dashed lines and each part represents a planar diagram. For each part, 
we define a certain quantity in order to show fl5.17p . For instance, if we consider the left 
diagram shown in FigJ21 we define the following quantity for the part labeled by s, 

A/2 A/2 A/2 A/2 ^ ^ ^ ^ 

5Z H H H (s-tfit-uY^iu-vYiu-wY' ^^'"^^ 

t=-k/2u=-K/2v=-K/2w=-K/2 ^ ' ^ ' ^ ' ^ ' 

(t^s) (u^t) (v^u) (w^u) 

Namely, for a dashed line connecting s and t, we assign ^^J-^^a; where 21 is the number 
of fields contained in the vertex and Qg is the summation of the product of them with 
respect to the variables other than s. Then, we can see that Qs has a finite value and 
bounded from above by a certain constant. This is shown as follows. First, note the 
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following identity, 

A/2 

E r^7)i7 = 2C(2/) - C(2/, A/2 + s + 1) - ({21, A/2 - s + 1), (B.2) 

t=-A/2 ^"^ ^ 

where ({z) and ({z, a) are the zeta function and the generalized zeta function, respectively. 
They are defined by, 

oo ^ oo ^ 

n=l n=0 ^ ^ 

Because C{2l,K/2 + s + 1) + C(2/, A/2 - s + 1) > 0, 

A/2 A/2 A/2 ^ ^ ^ 

«.< E E E (,_i)S(i_„)io(,_„)4 x2C(4) 

t=-A/2u=-A/2t;=-A/2 ^ ' ^ ' ^ ' 

(t^s) (n^t) (d^m) 



A/2 A/2 



t=-A/2 M=-A/2 



A/2 

< E 773^ 23c(4mio) 

t=-A/2 

< 24C(4)^C(8)C(10). (B.4) 

Thus, we find that Qs is finite even if we take the limit A — )■ oo. We are able to define 
Qs also for more general diagram shown in Figj3l If the diagram is 'tree', through the 
same calculation just as we described above, we see that Qs is bounded from above by a 
A- independent constant which is given by a product of the zeta functionq^. 

Then, let us consider the diagram in FigEl We can define Qs for this diagram. We 
assume that the rectangular part in Figj6] represents a 'tree' diagram so that Qs for this 
diagram is bounded from above. We consider the following quantity, 

A/2 A/2 A/2 

ATI E E 7737^7^^ = ATI E (C(2/)-C(2/,A/2 + s + l))Q.. (B.5) 

s=-A/2i=-A/2 ^ ' s=-K/2 



^Actually, the upper bound of Qs exists also for non-'tree' diagrams. However, this fact is not impor- 
tant for our purpose. 
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Figure 6: The rectangular box represents a Feynman diagram with 'tree'-hke structure, 
namely, a diagram which does not contain any loop formed by dashed lines. The circular 
blob represents a planar diagram which does not include dashed lines. 



Since there exists a positive A-independent constant c such that Qs < c, the second term 
on the right hand side in fIB.Sp is bounded as follows, 

A/2 A/2 

-— J2 C(2/,A/2 + s + l)g, <^ C(2/,A/2 + s + l) 

A+1 



A + _ , _ 

s=-A/2 s=-A/2 



2c 



A , _ 

n=l ?i=A+l 

where we have used (IB. 31) to obtain the second line. Because the second line in (IB. 6 1) is 
O (log A/ A) for Z = 1 and (9(1/ A) for / > 1, this quantity goes to zero in the limit A — )■ oo. 
Therefore, we find from ( ]B.5|) and (]B.6|) that 

lim -±- y y —^Qs = lim 4^ y Qs. (B.7) 

s=-A/2i=-A/2 ^ ^ s=-A/2 

Then, for instance, we can show the following equality by using f lB.7p iteratively, 

^ A/2 A/2 A/2 A/2 A/2 ^ ^ ^ ^ 

a^aTT 5Z 5Z 5Z (s-t)8(t-M)io (M-t;)4 (m-m;)4 

s=-A/2t=-A/2«=-A/2i>=-A/2i«=-A/2 ^ ^ ^ ^ ^ ^ ^ ^ 

(t^s) (v^u) {Wf^u) 

^ A/2 A/2 A/2 A/2 ^ ^ ^ 

"a^ATT Z Z Z Z (s-t)8(t-M)10 («-^;)4 ^ ^^^^^ 
s=-A/2t=-A/2M=-A/2i)=-A/2 ^ ' ^ ' ^ ' 

(t^s) (M^t) («^m) 

^ A/2 A/2 A/2 ^ ^ 

= £™ ATT E E E (,_y(,_„).o X2W 

s=-A/2t=-A/2M=-A/2 ^ ^ ^ ^ 
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A/2 

= i"^ A-T E 2^C(4)^C(8)C(10) 

A-5-00 A + i ^ — ^ 
s=-A/2 

= 2t(4)2C(8)C(10). (B.8) 

This is nothing but f l5.16p . Since (IB.SP imphes that we can replace Xlfi-A/2 ^/('^ ~ ^)^' 
on the tip of a branch in the 'tree' diagrams with C(20 the A — > oo hmit, applying the 
same calculation, we can show fl5.17p for the generic 'tree' planar diagram shown in FigJS] 



C Perturbative calculation 

In this appendix, we explicitly evaluate the free energy of (15 ■4p and the VEV of f l5.23p up 
to C(A'^). We perform the usual perturbation theory in (15 ■4p . 

First, we will calculate the free energy. Let us expand the interaction terms in terms 
of the power of the coupling 

oo 

^^('^^) - E 2Us-m 4^ E ^.C^^trCtrl-^,)^^-'", (CI) 
and the partition function as 

oo 

Z = Y,Z^. (C.2) 

fc=0 

where is the part of which the coupling dependence is 0[}^\ Zq is the free part, 

2„ = /n*«"^^-"*'=(|) . (C.3) 

and Zfc (fc = 1, 2, 3) are given by 

|i = -(V-.>. 

^0 

f = -(ly + \{v;'), 
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Figure 7: Feynman diagrams corresponding to (lC.7p . The left is planar diagram {0{N^)), 
while the right is nonplanar ((9(1)). 



§^=-{Vs) + {V2V,)-^{V,'). (C.4) 



We define the free energy of our matrix model as 

= In — 

2o 



ln(l + |l + ^ + |^ + ... ). (C.5) 



^0 ^0 ^0 

Substituting ( 1C.4I) into ( 1C.5I) . we obtain 

J" = -{Vl)c - {V2)c + ^(V^l\ - (V^3)c + (^2^l)c - ^(V'l')e + • " • , (C.6) 

where (■ ■ ■)c means the connected part of (■■■). For example, {Vi)c is calculated as follows; 



1 2, 1 / . 1 / A ^ ^ 



^ Q -f- ')! 



2 N 2 V ^ 2(s - t)2 ■ 2(t - m)2 V 47r2 



s,t,u 

1 / A \^ 

^'iV2(A+l)C(2)2 + ^(A+l)C(4) 



4vr4 ^ ^ Svr^ 

= (A + 1) AT^— + — . (C.7) 

In Feynman diagram introduced in subsection 15.11 one can describe this contribution as 
Fig. [3 Performing the similar calculations, we finally obtain f IS.lSp . 

Second, we perform a perturbative calculation of the VEV of the unknot Wilson loop 
fl5:23|l up to 0{}?). We expand the VEV as 
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1 °° u 

= 1 + (XTTW (20! + M<pru), (C.8) 

where (■ ■ ^ means the ladder part and (■ ■ ■)int means to evaluate with interaction vertices. 
The ladder part is easily calculated as follows 

J](tr0^)o = (A + i)iV, 

s 

s ^ ^ 

J2{tT<P% = {A + I) (5iV+^)- 

Then, up to C(A^) the ladder part gives 

1 >^ A' ^,^2. ^ A A2 / 1 \ A3 / 10 A 



(C.9) 



(A + 



(C.IO) 



Next, let us expand the third term in (IC.SI) up to 0{\^) 



1 °° A' 



mi 



(A + l)iV 



A 



(C.ll) 



For example, we can calculate the second term in (IC.lip as 



A= 



327r4 

SA^A^ 
' 47r4 



(2 X Ar(tr02tr0^)^ + 6(tr02tr02)^(tr02)^ x 2) 



1 



2(s -t)^ 



-(A + l)Arf 



1 + 



(A + 1)C(4) 



A^ 



The Feynman diagrams corresponding to the leading part are depicted in Fig. [HI We can 
calculate the other terms in the same way and the result is 



{A + 1)N {2iy. ^ ' 



int 
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Figure 8: Wilson loop corresponding to the leading part of (1C.12|) . 



A /A / 1 A A= A^ 1 A^ ^ / 1 A 2A 



Finally, gathering (IC.lOp and (IC.lSp . we obtain perturbative expansion of the unknot 



Wilson loop in our matrix model fl5.26p . 
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